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Summary. — We show that the excess of high energy neutrinos observed by 
the IceCube collaboration at energies above 100 TeV might originate from baryon 
number violating decays of heavy shadow baryons from mirror sector, which in turn 
constitute Dark Matter. Due to tiny mixing between mirror and ordinary neutrinos, 
it is possible to explain the specific features of the IceCube events spectrum. 

PACS 14.60.Pq - Neutrino mass and mixing. 

PACS 14.60.St - Non-standard-model neutrinos, right-handed neutrinos, etc.. 


Discovery of high-energy neutrino events by the IceCube Collaboration [1] opened a 
new era of experimental high-energy neutrino astrophysics. Recently the collaboration 
published the three year data collected between 2010 and 2013, containing 37 candidate 
in the energy range from 30 TeV to 2 PeV. While the number of events with energies 
below 100 TeV or so is compatible with the background expectations, there is an evident 
excess of the events with E > 100 TeV or so. On the other hand, no events were observed 
in the gap between 400 TeV and 1 PeV, but three most energetic shower events emerged 
at the end of the spectrum with energies 1.04 PeV, 1.14 PeV and 2.0 PeV where the 
atmospheric background is practically vanishing. On the other hand, the spectrum is 
apparently cut off at energies larger than about 2 PeV. The gap in the energy spectrum 
is difficult to explain in known models of high-energy neutrinos of astrophysical origin. 

The model that may explain such a spectrum was suggested in ref. [2]. Here we 
present main guidelines of this model. Nowadays it becomes popular the idea that 
dark matter of the universe can emerge from a parallel gauge sector, with particles and 
interactions sharing many similarities with ordinary particle sector. Such a shadow sector 
would contain particles like quarks which form composite baryons, as well as leptons and 
neutrinos which are all sterile for ordinary gauge interactions. In particularly interesting 
example is represented by so-called mirror world [3], which has the content of particles and 
interactions exactly identical to that of ordinary sector, with the same gauge and Yukawa 
coupling constants. Taking into consideration also attractive possibilities for particle 
physics beyond the Standard Model related to supersymmetric grand unified theory 
(SUSY GUT), one can consider that at higher energies both ordinary and mirror sectors 
are presented by identical GUTs, e.g. SU{5) and S'[/(5)' or SU{6) and S'C/(6)', which 
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in both sectors break down to their Standard Model subgroups SU{3) x SU{2) x 17(1) 
at the scale Mq ~ 10^® GeV. However, following refs. [4, 5], one can assume that 
later on the symmetry between two sectors is broken so that the electroweak symmetry 
breaking scale V in mirror sector is much larger than than in the Standard Model, 
V = 174 GeV. Namely, if V ~ 10^^ GeV, then the lightest shadow baryons should 
have mass order few PeV. Due to baryon violating GUT gauge bosons, they decay with 
Tdec ~ tu producing energetic shadow neutrinos which oscillate into active neutrinos 
(with oscillation probabilities ~ 10“^° or so) transferring their spectrum to the latter. 

Let us consider, for simplicity, a supersymmetric grand unification theory SU{5) x 
SU{5y. (^) Consider a situation when in both sectors the GUT symmetries are broken 
at the scale Mq ~ 2 x IQi® GeV, SU{5) SU{3) x SU{2) x U{1) and S'[/(5)' ^ 

SU{3y X SU{2y x U{iy. Below this scale our sector is represented by MSSM with chiral 
superfields of quarks qi = {u,d)i, u^dy and leptons li = {v,e)i, e^, i = 1,2,3 being 
family index, and two Higgs superfields hi,h 2 - As for parallel sector we have the similar 
particle content, quarks Qi = {U, D)i, Uf, Df, leptons Li = (TV, E)i, Ef, and two Higgses 
i7i,i72. At the scale /i = Mq all gauge coupling constants are equal and the Yukawa 
couplings have the same pattern for ordinary and shadow fermions: 


( 1 ) + Vy^q^d^hi + Yy^q,u^jh2) + (Yy^L.EyHi + Yf^Q.D^Hi + Y.yQ^UyH^) 

As for the neutrinos, their masses emerge from the following D = 5 operators [4]: 


( 2 ) ^ilh2)ilh2) + ^{LH2)iLH2) + ^{lh2){LH2) 


where M is an effective scale of the order of Planck scale Mpi [5]. The first two terms give 
the Majorana masses respectively to ordinary neutrinos and their shadow (sterile) 

partners while third term induces the mixing between ordinary and shadow neutrinos. 

Below GUT scale, the gauge coupling constants of two sectors, respectively ( 73 , 2,1 
and 53 2 1 evolve down in energies in both sectors in the same way (see Fig. 1). For 
ordinary sector we follow to the standard paradigm. Supersymmetry is broken at the 
scale Ms ~ 1 TeV, the Higgses are not protected anymore to be massless and they get 
VEVs Vi = vcos/3, V 2 = vsinfd, where v = 174 GeV, which induce the electroweak 
symmetry breaking and generate the fermion masses. 

As for parallel sector, let us assume now that mirror symmetry between two sectors 
is spontaneously broken in the following way. Imagine that supersymmetry is broken in 
shadow sector at the scale M'g ~ 10^^ GeV, supposedly due to non-zero E or D terms 
of some auxiliary fields, so that shadow scalars, including squarks and leptons as well 
as Higgs doublets i7i and 772 acquire soft masses order Ms-i^) Respectively shadow 
Higgses Get VEVs Ui = Ucos/3' and V 2 = Usin/3', and electroweak symmetry is broken 
at the scale V < M^. Therefore, the masses of shadow fermions are rescaled, modulo 
renormalization factors order 1, by a factor ^ = Vjv with respect to ordinary fermion 


(^) As discussed in ref. [2], our proposal can be more nicely realized in SUSY SU{6) theory [6] 
which gives natural solution to the so called hierarchy and doublet-triplet splitting problems by 
relating the electroweak symmetry breaking scale to the supersymmetry breaking scale Ms ~ 
1 TeV and naturally explains the fermion mass spectrum. 

(^) Interestingly, if supersymmetry breaking is transferred to our sector via gravity or other 
Planck scale mediators, this would explain ordinary soft masses order Ms ~ M'^/Mpi ~ 1 TeV. 
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Figure 1. - Running of gauge couplings below the GUT scale, Qi = 5i^/47r and a' = g'^/dvr. 


masses. Namely, taking V/v = 10® and assuming tan/3' = tan/3, by the RG running 
of gauge and Yukawa constants from the GUT scale down in energies, one obtains for 
the shadow electron mass Me — 0.4 PeV, Me — 1.1 PeV and Mu — 1-9 PeV, against 
TOe = 0.5 MeV, md — 5 MeV and ~ 3 MeV for ordinary electron and quarks. One 
the other hand, the RG evolution shows that the shadow QCD scale becomes A' ^ 100 
TeV (c.f. A ~ 200 MeV in ordinary QGD). Therefore, A' ^ Mu, Me and the shadow 
QGD looks like a rescaled version of our QGD without light quarks, but only with heavy 
quarks like c and b. Notice also that in shadow sector up quark U becomes heavier than 
the down quark D [2], and Mu,e/^' is of the same order as mb^c/^ in ordinary QGD. 

The lightest state would appear to be a shadow A“ baryon of spin 3/2, consisting of 
three down quarks D and having mass M/\ « 3Me = 3.3 PeV. All states, containing up 
quark U, will be unstable against weak decays, U —)■ DEN. As for mesons, the lightest 
pseudoscalar is shadow neutral pion tt® consisting of DD, with mass Mq « 2Me = 
2.2 PeV, while the lightest vector meson p^{DD) is slightly heavier than tt®. Charged 
Pion TT~ as well as p“-meson consisting of DU will have mass M ~ Mu + Me = 3 PeV, 
with p~ a bit heavier than 7 r“. All pseudoscalar and vector mesons have excited states 
with mass gap order A' between the levels, just like cc or bb states in ordinary QCD. 

Now we come to the role of baryon violation and proton decay which is fundamen¬ 
tal prediction of the GUTs. The heavy gauge bosons of SU{5) with baryon violating 
couplings between quarks and leptons induce the decay of the lightest ordinary baryons 
(proton, or neutron bound in nuclei), with lifetime Tp ~ AdQ{a\mp)~^ ~ 10®^ Gyr or so, 
where 05 is SU{5) coupling constant at the GUT scale Mq ^ 2 x 10^® GeV [7]. 

In the shadow sector, the similar couplings of GUT gauge bosons should destabilize the 
shadow baryon. However, taking into account that shadow A baryon is much heavier than 
the ordinary proton, M^/rUp ~ 10®, its lifetime must be about 30 orders of magnitude 
smaller than the proton lifetime. Hence we get ta ~ ^ 10 Gyr or so, with 

lifetime comparable to the age of the Universe tu = 14 Gyr. 

The principal decay mode of A baryon is in vector mesons, A“ —>■ p~ +Nx, where the 
neutrino state is generically a superposition of shadow neutrino flavor eigenstates. 
However, taking into account the active sterile mixing ~ j5v/aV due to operators (2), it 
contains also small admixture of ordinary neutrinos. Each decay produces monoenergetic 
neutrinos Vj, with energies Ea = ^M^{1 — where Mq is the mass of p~ meson 

and are the masses of its excitations. In Fig. 2(a) the spectrum of neutrinos 

produced by decay of galactic dark matter is shown by sharp peaks (solid blue) for 
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Figure 2. - (a) Shadow neutrino spectrum produced by A-barion decay and subsequent decay 
of shadow pions (b) Spectrum of VHE neutrino events as predicted in our model. 


Mo,Mi,M 2 respectively being 3.0, 3.1 and 3.2 PeV. Due to close degeneracy between 
the masses of A-baryon and p~ mesons, the neutrino energies are Ea 1 PeV while the 
most of initial energy = is taken away by mesons p~. 

Vector mesons promptly decay into mirror pion and photon, p~ —>■ 7 r“ + 7 ', and 
subsequent decay of the pion produces the neutrino once again (solid red curves in Fig. 
2(a)). Shadow 7 r“ has two decay modes, two body tt~ ENe and three body tt~ 
tt^ENe- Interestingly, their branching ratios are comparable which fact is intimately 
related to the value A' ^ 100 TeV [2]. Two body decay produces neutrinos with a 
narrow energy spectrum concentrated around M^/2 ~ 1.6 PeV, while the three body 
decay, due to smaller phase space, produces less energetic neutrinos with a wide spectrum 
extending up to the value i?max = — Mq — Me — Mu — Me — Me — 0.4 PeV. 

Fig. 2(a) shows the final spectrum including the neutrinos produced by the decay 
of dark matter in the galactic halo, and extragalactic neutrinos produced by the decay 
of cosmological dark matter at large redshifts. The fraction of extragalactic neutrinos 
strongly depends on the decay time ta (in Fig. 2 we take ta = lOtf/). Fig. 2(b) shows 
how such a spectrum will be seen by the IceCube provided that active-shadow neutrino 
mixing angles are order 10“^. Here the effective areas for the neutrino detection by 
IceCube [1] and characteristic error bars in estimation of neutrino energies (of about 
13%) are taken into account. Needless to say, that the obtained spectrum of events look 
very much like the spectrum observed by the IceCube [1]. The validity of our model will 
be tested with increasing statistics by the IceCube collaboration. 
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